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Spinning-black-hole binaries: The orbital hang up 
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We present the first fully-nonlinear numerical study of the dynamics of highly spinning black-hole 
binaries. We evolve binaries from quasicircular orbits (as inferred from Post-Newtonian theory), 
and find that the last stages of the orbital motion of black-hole binaries are profoundly affected by 
their individual spins. In order to cleanly display its effects, we consider two equal mass holes with 
individual spin parameters Sjm^ — 0.757, both aligned and anti-aligned with the orbital angular 
momentum (and compare with the spinless case), and with an initial orbital period of 125M. We 
find that the aligned case completes three orbits and merges significantly after the anti-aligned case, 
which completes less than one orbit. The total energy radiated for the former case is « 7% while 
for the latter it is only ~ 2%. The final Kerr hole remnants have rotation parameters ajM = 0.89 
and ajM = 0.44 respectively, showing the unlikeliness of creating a maximally rotating black hole 
out of the merger of two spinning holes. 

PACS numbers: 04.25.Dm, 04.25. Nx, 04.30.Db, 04.70.Bw 
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I. INTRODUCTION 

Spinning black holes play an important role in some of 
the most energetic astrophysical phenomena in the uni- 
verse. They form part of the main engine of gamma-ray 
bursts, being much more efficient at converting matter 
into radiation than non-spinning black holes. They are 
also responsible for the radio jets observed in active galac- 
tic nuclei, and the merger of two non-aligned spinning 
black holes is the likely explanation for the rapid direc- 
tional changes observed in these jets when galaxies col- 
lide ^] . Recent estimates 12] of the spin of stellar mass 
black holes by spectral fitting of the X-Ray continuum 
set the rotation parameter of two dynamically confirmed 
black holes at ajM ~ 0.75. Accretion, of course, can spin 
up black holes, reaching up to a sub-maximal spin rate 
of ajM ~ 0.95, when magneto-hydrodynamics is taken 
into account |3, |j| . Other models using the combined 
effects of gas accretion and binary-black-hole coalescence 
suggest that black holes may be rapidly rotating in all 
epochs [3. 

Recently new numerical techniques to solve the field 
equations of General Relativity have been developed jo, 
0, S S HOl that make it possible to stably evolve black- 
hole binaries for several orbits and to cornpute the cor- 
responding gravitational waveforms [III llJ; HjI' Nu- 
merical simulations of unequal-mass black-hole binaries, 
along with the calculation of the merger kicks, have 
been reported in Refs [ij) llM H^ ■ While research has 
been mainly focused on initially non-spinning black holes, 
there are important questions to be addressed when we 
consider highly-spinning black holes (see J17l | and refer- 
ences therein). 

In this paper we study how the emission of gravita- 
tional radiation affects the orbital trajectory of highly- 
spinning, equal-mass black holes as a function of the 
spin orientation. In order to maximize the effect, we 



consider black-hole binaries with both spins aligned and 
anti-aligned with the orbital angular momentum, as well 
as the corresponding spinless case. We shall consider 
quasicircular orbit initial data with the same initial or- 
bital period (as determined by the third post-Newtonian 
(3PN) expansion). In this way differences in the subse- 
quent evolution can be attributed to the differences in 
the generation and emission of gravitational radiation. 

In Ref. [l8J the numerical evolutions of spinning bina- 
ries were studied for relatively modest values of the spins 
(—0.25 < S/m? < 0.17, m being the horizon mass of the 
individual holes) using the 'Lazarus' technique of match- 
ing full numerical evolutions to perturbation theory. In 
those evolutions the spin of the remnant Kerr hole in- 
creased with S/M"^ for the aligned case. Extrapolation 
to maximally spinning individual holes indicated that the 
remnant would remain sub-maximal for Sjm? < 0.85. 
We will revisit this scenario, now reaching much higher 
values of the individual spins in order to make a more 
accurate statement. 



II. INITIAL DATA 

We use the Brandt-Briigmann puncture ap proach 
along with the elliptic solver BAM_Elliptic [l a l20 | to 
compute initial data. Table ^ gives our choice of ini- 
tial parameters. We have taken a fiducial angular fre- 
quency of Mfl = 0.05, which corresponds to an orbital 
period of approximately T = 125M. This, accordingl y to 
our previous simulation for non-spinning black holes |ll| 
makes the binary complete more than a full orbit be- 
fore the black holes merge. We choose individual spins 
S — ±0.757 m^ (as measured using isolated horizon tech- 
niques |2l| ) to guarantee that the total angular momen- 
tum in the aligned case exceeds M^, the maximum al- 
lowed value for a common horizon to form. The gravita- 



TABLE I: Initial data for quasicircular orbits of black-hole bi- 
naries with spin. The holes have proper horizon separation I, 
with puncture locations (0, ±y, 0), linear momenta {^P, 0, 0), 
and spin (0, 0, 5). J is the total angular momentum, L is the 
orbital angular momentum, f2 is the orbital frequency, nip is 
the individual puncture mass. All in units of the ADM mass 
M. 



S/m' 


-H-hO.757 


0.00 


-- 0.757 


l/M 


9.27 


9.62 


10.34 


y/M 


3.0595 


3.280 


3.465 


P/M 


0.1291 


0.1336 


0.1382 


S/M^ 


+0.1939 


0.000 


-0.1924 


J/M"" 


1.1778 


0.8764 


0.5729 


L/M'' 


0.7900 


0.8764 


0.9577 


Mn 


0.0500 


0.0500 


0.0500 


nip/M 


0.3344 


0.4851 


0.3344 



tional radiation emitted should efficiently carry out angu- 
lar momentum from the system in order for the cosmic 
censorship conjecture to hold 22]. We can thus begin 
to explore its validity here, and this will be the subject 
of a more detailed study in a forthcoming paper by the 
authors. 

With our choices of the spins and the orbital angular 
frequency, we determine the remaining orbital parame- 
ters by imposing quasicircular orbits according to the sec- 
ond post-Newtonian expansion of spinning particles |23l | 
extended by the third-order orbital corrections 12411. We 
then give those parameters to the exact Bowen J23 so- 
lution of the momentum constraints and solve for the 
conformal factor of the (conformally flat) three-metric 
to complete our choice of the initial data. This post- 
Newtonian data should produce orbits with acceptably 
small eccentricities, as can be seen when comparing the 
zero-spin parameters in table ^ with others proposed in 
the literature (e.g. [2g). We have also briefly studied the 
effects of a different choice of the form of the initial data 
(Kerr conformal extrinsic curvature) for spinning black 
holes, as proposed in Ref. 27]. However, the spurious ra- 
diation in the initial data is dominated by the momentum 
terms and both data sets give nearly identical waveforms 
(see Fig. O. 



III. TECHNIQUES 

We evolved these black-hole-binary data sets using the 
LazEv 28] implementation of the moving puncture ap- 
proach 9. lOj . In our version of the moving puncture 
approach Q we replace the BSSN [23, Hfl Hj confor- 
mal exponent 0, which is infinite on the punctures, with 
the initially C* field x — cxp(— 40). This new variable, 
along with the other BSSN variables, will remain finite 
provided that one uses a suitable choice for the gauge. 

We obtained accurate, convergent waveforms by evolv- 
ing this system in conjunction with a modified 1+log 
lapse, a modified Gamma-driver shift condition |9l. l32| . 



and an initial lapse a ^ i^BL- ^^'^ lapse and shift are 
evolved with {dt - (3'd,)a = -2aK, dtP"" = B", and 
dtB'^ = 3/AdtT°- — rjE"^. These gauge conditions re- 
quire careful treatment of x near the puncture in or- 
der for the system to remain stable |3, [ill- For our 
version of the moving puncture approach, we find that 
the product aA^^djCJ) has to be initially C* on the punc- 
ture. In the spinning case, A*-' is 0{r^) on the punc- 
ture, thus requiring that a oc r"^ to maintain differen- 
tiability. We therefore choose an initial lapse a ~ "(Pbl 
(which is 0{r'^) and C^ on the puncture). In particu- 
lar, a{t = 0) = 2/(1 + V'bl) reproduces the isotropic 
Schwarzschild lapse at large distances from the hole. The 
initial values of /3* and B^ were set to zero. 

The minimum resolution required to accurately model 
the dynamics of the merger scales with rup. We would 
expect satisfactory results for a minimum resolution of 
h — M/30 (based on the non-spinning case, where sat- 
isfactory results were obtained with h = M/21, and the 
ratio of the puncture masses in the spinning and non- 
spinning cases), M being the total ADM mass. However, 
the additional power of 1/r introduced to Kij because 
of spin, necessitates even higher resolution (we estimate 
M/40 — M/50) to get highly accurate waveforms. 

We exploited the Pi-rotational symmetry about the z- 
axis as well as the refiection symmetry about the xy plane 
to reduce the size of the computational domain by one 
quarter. 



IV. RESULTS 

We evolved the '--0.757' configuration using grid sizes 
of 320^ X 640, 3842 ^ jQg^ ^^^ 4432 ^ ggg and resolutions 
of Af/25, M/30, and M/35 respectively. We used a mul- 
tiple transition fisheye [Tl| to push the physical bound- 
aries to 134M. We calculate t/ji in the Quasi-Kinner sley 
frame using the recently developed techniques of Ref. |33| 
that allow a meaningful extraction closer to the hole. In 
Fig. n we show the real part of the {£ = 2,m = 2) mode 
of rip4 for the -- 0.757 case (extracted at r = lOM) for 
these resolutions, as well as a convergence plot of these 
data. The waveforms show fourth-order convergence up 
to i ~ llOM. The phase error from the h = M/25 run 
becomes too large to measure a meaningful convergence 
rate after t ~ llOM. Higher resolution runs will remain 
convergent, as demonstrated by the better phase agree- 
ment between the M/30 and M/35 runs. We extract the 
waveform at lOM to minimize the effects of the extreme 
fisheye deresolution (which is too strong in the far field 
to get accurate waveforms). After a time translation, 
the waveforms from the three runs only differ by a con- 
stant phase factor. We calculate this factor and plot the 
phase-corrected waveforms [la Isi) in the upper inset of 
Fig. n Note the near perfect agreement both in the or- 
bital and plunge waveforms. The waveforms calculated 
at r = lOM do not yield accurate estimates for the ra- 
diated energy (as expected, since the observer is still in 
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FIG. 1: The real {£ — 2,m = 2) component of ri/)4 in the 
Quasi-Kinnersley tetrad at r = lOM for the -- 0.757 case. 
The lower inset shows the differences rtp4,{M /25) — ril)4{M /30) 
(solid line) and ri/)4(M/30)-rV'4(M/35) (dotted line), the lat- 
ter rescaled by 2.33 to demonstrate fourth-order convergence. 
The lack of convergence for t < lOM is due to roundoff effects 
in the initial data solver. The upper inset shows the real part 
of the phase-corrected {£ = 2,171 — 2) mode of '(/14 at the same 
radius. Note the near-perfect agreement after t = 45M. 



the near zone). However, as shovifn below, we obtained 
highly accurate estimates for the radiation be examin- 
ing the remnant horizon. In order to obtain accurate 
measurements for the radiated energy and angular mo- 
mentum from the waveform, one needs to use a weaker 
fisheye deresolution in the outer region, and carefully ad- 
just the gauge so that the waveform is highly accurate at 
large radii (r ^ 50M). 

The relatively large phase errors in this spinning case 
compared with our zero-spin simulations .dl | are due to 
the fact that the effective resolution in the spinning case 
is smaller due to the smaller value of nip as well as the 
lower order differentiability of the spinning data com- 
pared to zero-spin. A likely explanation is that numer- 
ical dissipation (which more strongly affects this higher 
frequency data) causes the merger to happen sooner than 
expected. This dissipation decreases with resolution. 

We used Jonathan Thornburg's AHFindcrDirect 
thorn ^3 to calculate the apparent horizons. We find 
that the common horizon is first detected at t = 105.5A/ 
and has a mass of M-h = 0.978 ± .OOIM and rotation 
parameter of a/Mu = 0.443 ± 0.001. During the merger 
(2.2 ± 0.1)% of the mass and (26 ± 2)% of the angular 
momentum are converted into radiation. 

Mu and ajM-yi were obtained from the asymptotic val- 
ues of the horizon surface area and the ratio of the polar 
to equatorial circumferences (see Refs [1,|33,|3^)- The 
ranges given for these quantities arise from the uncer- 
tainties in obtaining these values at finite time, and are 




FIG. 2: The puncture trajectories on the xy plane for the 
'--' case with resolution M/35. The spirals are the puncture 
trajectories with ticks every lOM of evolution. The dot-dash 
'peanut shaped' figure is the first detected common horizon at 
105. 5M. The (extrapolated) period of the last orbit is around 
120M. 



independent of resolution. Thus, these horizon parame- 
ters give an accurate and robust measurement (even in 
the unresolved h = M/2b case) for the radiated mass and 
angular momentum. 

FigureElshows the trajectories of the punctures for the 
- - 0.757 configuration as well as the projection of the first 
common horizon on the xy plane. It is evident from the 
waveform and the track that the binary undergoes ~ 0.9 
orbits before merging. 

We evolved the '-I--I-' configuration with a grid size of 
384^ X 768 and resolution of M/30. We used multiple 
transition fisheye to push the boundaries to 159M. In 
Fig. Owe show the (£ = 2,m = 2) mode of r'i/'4 in the 
Quasi-Kinnersley frame for the '-I--I-0.757' case, again ex- 
tracted at r = lOM. Note the 'plunge' waveform is de- 
layed by ~ 120M compared to the '--0.757' case. The 
waveform shows approximately six periods of orbital ra- 
diation prior to the plunge waveform, indicating that the 
binary completed approximately three orbits. 

We repeated the '-I— fO.757' case with a gridsize of 
448^ X 896 and resolution of M/30 to force the bound- 
aries to 266M. This new configuration allows us to ac- 
curately obtain the horizon parameters (since they are 
not contaminated by the boundary), but is too coarse in 
the far-field region to produce accurate waveforms. The 
first common horizon was detected at t = 232. 5A/. In 
this case the final horizon had a mass of 0.933 ± .OOlAf 
and spin of 0.890 ± .002 (indicating that (6.7 ± 0.2%) of 
the mass and (34 ± 1)% of the angular momentum are 
radiated away). Table HTl gives a summary of these re- 
sults. Note that the above values for the radiated energy 
are in rough agreement with those estimated using the 
effective one body approximation for maximally spinning 
holes 113. 

Figure 01 shows the track for the '-I-+0.757' configura- 
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FIG. 3: The real part of the {£ = 2,m = 2) mode oi rtpi in 
the Quasi-Kinnersley frame at r = lOM from the '+ + 0.757' 
configuration. (The small circles are the early-time waveform 
from conformal Kerr data.) The top inset shows a magnified 
view of the early orbital motion. Note that the '++0.757' 
waveform has 6 wavelengths of orbital motion prior to the 
plunge waveform (at t ~ 232. 5M), indicating that the bi- 
nary orbited approximately three times before merging. The 
bottom inset shows the real (solid) and imaginary (dotted) 
components of the (2,2) component of the strain h calculated 
at r = lOM. 



tion. Note that the spiral is much tighter than in the 
'--0.757' configuration, and that the binary completes 
roughly 3 orbits before the common horizon forms. Note 
also that the first common horizon is much smaller in this 
case (in these coordinates) than in the '--0.757' case. 

To demonstrate consistency with the General Relativ- 
ity field equations, we calculated the Hamiltonian con- 
strain violation. The constraint converges to fourth-order 
outside a small region surrounding the puncture (the con- 
straint violation on the nearest neighboring points to the 
puncture is roughly independent of resolution, but this 
non-converging error does not propagate outside the in- 
dividual horizons). Figure |S1 shows the Hamiltonian con- 
straint violations for the '--' configuration along the x- 
axes at t = 45M and along the y-axis at t = 80M (at the 
time when the punctures cross the x-axis and 5M after 
the punctures cross the y-axis for the second time) for 
the M/30 and Af/35 runs. The constraint is convergent 
everywhere except points contaminated by boundary er- 
rors (these points have been removed from the plot) and 
at the points closest to the puncture. 

We complete our initial study with the corresponding 
spinless case as a reference point. For details on the accu- 
racy and evolution of the spinless case see Ref. [ij . We 
evolved the zero-spin case with a resolution of M/22.5 
and gridsize of 320^ x 640 (the outer boundary was lo- 
cated at 216A/). The first common horizon formed at 
t = 161M with mass (0.965 ± .001)% and spin a/M = 
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FIG. 4: The puncture trajectories on the xy plane for '-I— 1-' 
configuration with resolution A//30. The spirals are the punc- 
ture trajectories with ticks every lOM of evolution. The dot- 
dash 'peanut shaped' figure is the first detected common hori- 
zon at i = 232. 5M. The period of the last orbit is around 
36Af . The last orbit begins when the punctures are located 
at 1.4Af from the origin (in these coordinates). 




FIG. 5: The Hamiltonian constraint violation ai t — 45 Af 
along the x-axis (top plot) and at i = 80Af along the y- 
axis (bottom plot) for the Af/30 and Af/35 runs (the latter 
rescaled by (35/30)*) for the '--' configuration. The punc- 
tures crossed the x-axis at i = 45M and crossed the y-axis 
for the second time at i = 75M. Note the reasonable fourth- 
order convergence (except at the puncture). Points contami- 
nated by boundary errors have been excluded from the plot. 
The high frequency violations near the numerical coordinate 
y/M = ±9 are due to the extreme fisheye deresolution near 
the boundary, and converge with resolution. 



(0.688 ± .001). This corresponds to a radiated energy and 
angular momentum of (3.5 ± 0.1)% and (26.9 ± 0.1)% re- 
spectively. 



TABLE II: Results of the evolution as determined from the 
remnant horizons. The horizon formed at f = 224. 5M for 
the '++0.757' configuration, t = 105. 5M for the '--0.757' 
configuration, and t — 161A/ for the spinless configuration. 



S/r 



Erad/MADM Jrad/ JaDM 



a/M', 



H 



-H+0.757 

-- 0.757 

0.00 



(6.7 ±0.2)% 
(2.2 ±0.1)% 
(3.5±.1)% 



(34 ± 1)% 
(26 ± 2)% 
(26.9 ±0.1) 



0.890 ± 0.002 
0.443 ± 0.001 
0.688 ± 0.001 



V. DISCUSSION 

In this paper we have shown that the 'moving punc- 
ture' approach can be used to accurately simulate the 
inspiral orbit of spinning-black-hole binaries. We found 
that the spin-orbit coupling delays the onset of the plunge 
phase (compared to the non-spinning case) when the 
spins are aligned with the orbital angular momentum, 
while in the anti-aligned case the plunge phase is has- 
tened. In all cases, the black holes merge to form a single 
Kerr black hole with rotation parameter a/M < 1. 

A fit of the remnant spin to the spins of the initial 
black holes leads to the simple extrapolation formula 
a/M « 0.688 + 0.298{S/m^) - 0.038{S/mY- Although 
more accurate simulations are needed, these results show 
that it is very unlikely to form a nearly maximally rotat- 
ing black hole out of the merger of two highly spinning 
ones. Our results reinforce the same qualitative conclu- 
sions reached with the Lazarus approach [l3 |. and are 
consistent with those in Refs [l7l ISTJ . Extrapolation of 
the radiated energy to maximally rotating black holes 
with the fit Erad/M « 3.5 + 3{S/m^) + 152/90{S/m^ 
set it near 8.2%, not far from the assumed 10% in Ref. [13 
a decade ago. 

For the '--0.757' case the final orbit lasts ^ 120A/ 
starting at a separation of 7M in coordinate space, while 
for the '++0.757' case we have found that the duration 
of the last orbit is ~ 36M at a coordinate separation 
of 2.8M. It is worth noting from "Sg, that the orbital 
period associated with the ISCO for the '++0.17' case 



is roughly 27M at a coordinate separation of 1.6M, and 
that no ISCO was found for higher spin configurations. 
This highlights again the importance of the gravitational 
radiation in the late binary black hole dynamics which is 
not captured in the determinations of the ISCO. On the 
other hand, the dependence of the ISCO on spin correctly 
implies that the '++' configurations are more stable at 
close separations than the '- -' configurations. This stabil- 
ity property is observed in the significantly tighter spiral 
displayed in Fig. ^ 

The Post-Newtonian equations of motion |23| indicate 
that the leading spin-orbit interaction is of 1.5PN order, 
while the spin-spin interaction is of 2PN order. It is 
the spin-orbit interaction (attractive/repulsive for --/±+ 
configurations respectively) responsible for the longer the 
stability of the aligned spin binary. 

Many outstanding issues involving spinning black hole 
scenarios remain to be explored. We plan to study some 
of them next, including additional values of the individ- 
ual spins for the '++' cases in order to better extrapolate 
the results to the limiting maximally rotating individual 
holes, as well as unaligned spins to study precessional ef- 
fects. Finally, more significant computer resources and 
the use of Adaptive Mesh Refinement techniques com- 
bined with higher order finite difference methods will be 
needed to achieve the accuracy required to aid gravita- 
tional wave detection efforts l39l. 
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